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Abstract. We present a numerical study of the self-similar solutions of the Localized Induction 
Approximation of a vortex filament. These self-similar solutions, which constitute a one-parameter 
family, develop a singularity at finite time. We study a number of boundary conditions that allow 
us reproduce the mechanism of singularity formation. Some related questions are also considered. 
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1. Introduction. Given a curve Xq : M — > M'', consider the geometric flow 

Xt = cb, (1.1) 
where c is the curvature and b the binomial component of the Frenet-Serret formulae 



(1.2) 



The flow can be expressed as 
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(1.3) 



where A+ is the usual cross-product, and s denotes arclength. This flow appeared 
for the first time in 1906 ^lOj and was rederived in 1965 by Arms and Hama [T] as 
an approximation of the dynamics of a vortex filament under the Euler equations. 
This model is usually known as the Localized Induction Approximation (LIA). We 
refer the reader to [5 and [26, for an analysis and discussion about the limitations 
of this model. Starting with the work of Schwartz in [57], the LIA has been also 
used as an approximation of the quantum vortex motion in superfluid Helium. Of 
particular relevance for our purposes is the recent work of T. Lipniacki [24) . (25] . A 
rather complete list of references about the use of LIA in this setting can be found in 
these two papers. 



Some of the explicit solutions of (1.3 1 are the line, circle and helix. It is easy 



to see that the tangent vector T = X^ remains with constant length, so that we 



can assume that it takes values on the unit sphere. Differentiating (1.3 1, we get the 
following equation for T: 

T, =TA+T,,. (1.4) 
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This equation, known as the Schrodinger map equation on the sphere, is a particular 
case of the Landau-Lifshitz equation for ferromagnetism |23j and can be rewritten in 
a more geometric way as 

Tt=JD,T„ (1.5) 

where D is the covariant derivative and J is the complex structure of the sphere. 
Written in this way, ( 1 .4 1 admits an immediate generalization and we can change 
both its definition domain (considering for instance more variables) and its image 
(considering other more complex varieties). We will insist on the second possibility, 
choosing also the hyperbolic plane as the target space; in that case, the equation 
for T is 



TA_T,, (1.6) 



and, equivalently, for X, 



Xt = X,A_X,„ (1.7) 

with A_ defined as 

a A- b = (0263 - a3&2, 03^1 - ai&3, -(0162 - a2&i))- 
In this article, we study numerically the self-similar solutions of 

Xt = X, A± X,,. (1.8) 

where A± has been defined as 

a A± b = (0263 - 0362, 0361 - 0163, ±(0162 - 02^1)). 

Equivalently, a generalized version of the scalar product, denoted by o_|_, is given as 

a o_|_ b = aibi + 0262 i 0363. 

By using ±, we can consider simultaneously the Euclidean case, corresponding to + 
and the Hyperbolic case, corresponding to — . Equivalently, whenever we use =F, the 
— sign will refer to the Euclidean case and the + sign to the Hyperbolic one. Using 
the A± notation, the equation for X^ = T is now 

Tt = TA±T,,. (1.9) 

If T e H^, it is still possible to give a generalized version of the Frenet-Serret trihedron 
(1.2) for each point of the curve X, formed by T and two other vectors ei and 62- 
Indeed, a few calculations show that, for both cases, all the possible generalizations 
of ( 1.2 1 have the form 

a f3\ 

Ta (5 I • I ei I , (1.10) 

-s 




where T o± T = ±1, ei oj- ei = 62 oj. 62 = 1, T oj- ei = T oj- 62 = ei oj- 62 = 0; with 
regard to A±, we have T A± ei = 62, ei A± 62 = — T, ei A± 62 = ±T. Without loss 
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of generality, we can choose one of the coefficients a, f3 or 6 to equal zero. If we make 
(3 — and denote a = c and 6 = t, ( 1.10 ) becomes 



(1.11) 



In the Euclidean case, ei and 62 refer to n and b respectively. In the hyperbolic case, 
we can refer to c and r as the generalized curvature and torsion. In both cases, it is 
possible to recover X, except for some rigid motion, if c and r are known. 

The self-similar solutions of (1.8 1 have been studied in in the Euclidean case 
and in |TT] in the hyperbolic one -see also |2B for some related work. They are 
such that if X(s,t) solves (I.8I, so does A~^X(As, A^t). Therefore, taking A = t^^/^ 



and defining G(s) = X(s, 1), the self-similar solutions will be of the form 
X(s, t) ^ ti/2x(t-i/2s^ 1) ^ y^Gis/Vi). 



(1.12) 



Bearing this in mind, after some straightforward but tedious computations |19l 111], 
the family of solutions in which we are interested is defined by 



where cq is the family parameter and G' = T(s, 1) is the solution of 

Co 0^ 





(1.13) 



(1.14) 



Xs , ei and 62 being the components of the generalized Frenet-Serret formulae; 



the initial conditions for (1.141 are 



'G(0)-2co(0,l,0), 
T(0,1)- (0,0,1), 
ei(0,l) = (1,0,0), 

,62(0,1) = (0,1,0). 



(1.15) 



Finally, for cq = we define 



Xo(s,i) = s(0,0,l). 
It can be easily shown that, for an arbitrary t > 0, we have 



(1.16) 



(1.17) 



The one-parameter family of self-similar functions we have just defined satisfies the 
following theorem: 

Theorem 1.1. Given cq > 0, the Xc„ defined by (1.131, (1.15) and (I.I6I is a 
C°° solution of ([TS]), Vt > 0. 

Moreover, there are A^{co), A^(co), B^(co), B^(co) and a constant C such that 
(i) |Xco(s,i) - Ais(co)x[o,+oo)(s) ~ A."^ s{c4))x(-oo.o]{s)\ < C^/t. 
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(ii) We have the following asymptotics: 



G(s) = A^ (co) (^s ± 2^) - 4co^ + 0{1 



±00: 



62 



T(s) = A^(co)-2co- + 0(l/s^), 
s 

(ei - le^) = B^(co)e^^'/4e±'^o'°s|s| _^ oil/s), 



±00: 



±00; 



where the superindex j = 1 when s 00 and j = 2 when s — > —00. 
(in) = {A{,A{,Ai) and = {B{,B{,BI) are vectors, with o± A^ = ±1 
and 



Bl = Bl Bl = Bl, Bl = -Bl A^o±B-'=0, 

We have combined in this theorem the results for both the Euchdean and the Hyper- 
bohc cases, proved in [19] and [TT], respectively. Naturally, A^ , B-' and C are different 
in each case. 



It follows from this theorem that there exists a solution to (1.8) such that 
X(s,0) = A^sx[o,+oo)(s) + A^sx(_oo,o](s), 



(1.18) 



and, correspondingly for (1.9), 

T(s, 0) = Aix[o,+co)(s) + K^X{-oo.»]{s). 



(1.19) 



Observe that both ( 1.8 1 and ( 1.9 1 are time reversible, because if X(s, t) and T(s, t) are 



their respective solutions, so are X(— s, —t) and — T(— s, ~t). Therefore, we have two 





Fig. 1.1. Comparison between the theoretical evolution of'K{s,t) and a real experiment 
of a colored fluid traversing a triangular wing. The experimental result on the right has been 
taken from f^/. 

one-parameter families of regular solutions that develop singularities at finite time. In 
the case of X, we have precisely a corner-shaped singularity, having thus a good model 
to describe some natural phenomena. In Figure |1.1[ for instance, we have plotted on 
the left hand-side the evolution of X(s,i), obtained by integrating (1.17), and on the 



right hand side, a picture of an experiment published in ONERA [29^, where several 
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lines of colored fluid in water show the symmetrical pair of vortices behind an inclined 
delta wing, with a high Reynolds number. The resemblance between both images is 
really striking, at least at the qualitative level. In fact, in both images we clearly 
see the self-similar behavior of the evolution of the filaments. Also in both cases the 
filaments tend asymptotically to two non-parallel lines. Finally, both pictures have in 
common the oscillatory behavior, and the shift of the horse-shoe- type of curve that is 
close to the vertex with respect to the plane that contains the two asymptotic lines. 
See also figures 7.5.7 and 7.8.6 of plate 2 in ^Sj. It would be interesting to know if this 
resemblance is not just at the qualitative level. Notice that Theorem 1 1 . 1 1 quantifies in 
a precise manner the above mentioned properties for the self-similar solutions of the 
LI A. 

The main difference between the Euclidean and the hyperbolic case is the fact 

that 

A\ = Al^e^^^^ , (1.20) 
i.e., grow exponentially with cq in the hyperbolic case. Thus, although T, ei and 
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Fig. 1.2. Co = 2.7, 



olic case. Euclidean lengths of T, ei and 62 



62 are also bounded in the hyperbolic case for a given cq, this is by no means trivial, 
unlike in the Euclidean case; moreover, there are no global bounds valid for all cq and 
the bounds grow very fast even for relatively small values of cq . In Figure |1.2| for 
instance, we have plotted the Euclidean lengths of T, ei and 62 in the hyperbolic case, 
for Cq = 2.7; the lengths are seen to be 0(10^). This fact makes the proof of Theorem 
|l.l|m uch more difficult in the Hyperbolic well as the numerical treatment of 



(1.8) and (1.9) for large values of cq. 



Now, (1.8 1 is invariant under translations and also under rotations in the Eu- 



clidean case and under Lorentz transformations with unitarian determinant in the 
hyperbolic case; hence, given two arbitrary in S'^ or H^, we can transform (I.8I, 
so that A^ have the form as in part (iii) of Theorem |l.l[ Then, it is possible to 
obtain the corresponding cq. Observe that for all A-' G there is a corresponding 
Cq. However, when A^ = A^ G which can be reduced to A^ = A^ = (0,0, 1), Cq 
should be infinity, so there is no real co matching that case. 



There is a natural connection between equations (1.8 1 and (1.9), and the non- 



linear Schrodinger equation. Equation (1.8 1 can be transformed into the Schrodinger 



equation with a cubic nonlinearity via a transformation introduced by Hasimoto 
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Especifically, define 



ip{s,t) = c{s,t) exp [ i I T{s',t')ds' 



(1.21) 



Then, ij) satisfies the equation 



(1.22) 



for sonic function A{t), which can be removed by means of an integrating factor. 
Equation (1.9 1 can be transformed into a nonlinear Schrodinger equation by per- 



forming a stereographic projection onto the complex plane. If we define 



z{s, t) — , , + i ^ 



1 + T. 



then z satisfies the equation 



Zt = IZss T 



2iz 



l±\z 



zl 



2 s 



(1.23) 



(1.24) 



In this article, we study the self-similar solutions described above from a numerical 
point of view, in order to understand the singularity formation, as well as the mecha- 
nism for energy concentration that is responsible for the formation of such singularity. 
These solutions can be studied at the level of X (eq 
(1.22l), or at the level of z 



(eq. (1.241) 



(1.80, T (eq. (1.9)), V (eq 



(1.8), given that X can be recovered from 
is fixed by considering ( 1.13| ), i.e.. 



However, we will not work directly with 
^ except for a constant of integration that 



X(s,0) =2co\/t(0, 1,0). 



(1.25) 



Although there is a rich literature concerning the numerical study of ( |1.4| and, in 
general, the Landau-Lifshitz equation (see, for instance, O [HI [171 [13] ) , the materials 
analyzing the numerics of the self-similar solutions we are considering are very scarce; 
the most relevant results were given by Buttke some twenty years ago [7] . Nevertheless, 
Buttke only considered the forward case, starting with a singular initial datum at t = 
and going forward in time; we are interested in the singularity formation process, and 
therefore will focus mostly in the backward case, i.e., starting from t = 1, we will 
approach the singularity time. 

Moreover, we are not aware of any other study of the hyperbolic case. For small 
Co, the results of both cases are virtually identical; because of that, we have used in 
our experiments cq = 0.2. In contrast, for bigger cq, because of the rapid growth of 



the Euclidean lengths of T, ei and 62 (see Figure 1.2), the hyperbolic case becomes 



much more difficult to treat numerically, requiring a finer study, which we postpone 
for the future. 



The three equations (1.9), (1.22), and (1.24) are posed in the whole real line. To 
carry out the numerical simulations, we must restrict ourselves to an interval [— L,I/] 
(for L large enough) . It is important then to give appropriate boundary conditions in 
order to capture the singularity formation mechanism. In [T, only periodic boundary 
conditions were considered. 

The experiments performed here constitute numerical evidence of the stability of 



(1.8), (1.9) and (1.24), and in particular of the robustness of the energy concentration 
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mechanism. Indeed, one can never introduce numerically completely exact boundary 
conditions at s = ±L; in the best of the cases, there will be small errors due to machine 
precision and in the worst one, we will be using very rough boundary conditions. 
Therefore, we can consider that we are computing perturbed versions of the exact 
solutions. In a recent series of articles [3l H] Banica and Vega prove the stability 
of the solutions that we consider at the level of T. At the level of ip, Banica and 
Vega |4| also prove that a logarithmic instability appears when t comes close to zero. 
Notice that T involves one integral with respect to the curvature and two with respect 
to the torsion. In our case these are highly oscillatory integrals so that, due to the 
cancelations, T and X become stable. These results give the theoretical basis to the 
numerical results presented here, and suggest that it is better to work numerically 



with the equations at the level of T or z. Therefore, we consider only equations (1.9 1 



and (1.241 



We return now to the issue of boundary conditions. For equation (1.9 1, we con- 
sider two types of boundary conditions, derived from the asymptotics of T given in 
Theorem 1 1.1[ The first boundary condition is simply T{±L,t) = T(±, 1) and for the 
second one, we keep the first-order term in the asymptotics of T(s, 1) and use the fact 
that T{s,t) = T{s/Vi, 1). 

For equation ( |1.24[ ), we consider three types of boundary conditions: the first one 
is the projected boundary condition used for T, The second boundary condition is 
derived from the self-similarity condition, i.e., z{s,t) = z{s/y/i, 1). Unlike the former 
one, this condition is not an approximate one, so it introduces no noise, giving very 
clean results. Finally, we derive a radiation boundary condition. This condition seems 
to perform best, and in contrast to the previous two conditions, can be used in the 
progressive case. 



The remainder of the article is organized as follows: Equation (1.9 1 is considered 
in section [2] where we use a simple spatial discretization using finite differences. The 
nonlinearity in the equation and the constraint T=l, make it difficult to construct 
efficient implicit methods for these equations |16| . For the two boundary conditions 
mentioned earlier, we observe an stability constraint of \At\ = 0{As^). The first 
condition only reproduces the singularity formation from a qualitative point of view, 
while the second one gives good results for all s, even for very small t, although a very 
small As may be needed. 



In section [3] we consider equation (1.241. This equation adapts well to a pseudo- 



spectral method using Chebyshev polynomials; comprehensive information about this 
family of polynomials can be found, for instance, in [14 and [TB]. We use also a 
Chebyshev point distribution, which is now much better suited, since a denser con- 
centration of points is required near the boundary for small t and, specially, for big 
L, which is when the pseudo-spectral method really excels. 

In order to measure the accuracy of our results, we compute the error in the 
curvature obtained from T (which, when necessary, can be recovered from z as well), 
since the correct curvature has the known value c(s,t) — 

Finally, the stability results of Banica and Vega are described in section |4] 
2. Schrodinger map. We consider equation 

f Tt(s, t) = T(s, t) A± T,,(s, t), s G [-L, +L] 
1t(s,1) = T0(s). 
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In his Ph. D. Thesis, Buttke considered the following Crank-Nicholson type numerical 



scheme to integrate the Euclidean version of (2.1) 



T(s, t + At)- T(s, t) _ T(s, t + At)+ T(s, t) ^ f D+^T{s, t + At) + D+^T{s, t) 



At 2 V 2 

(2.2) 

where Dj, represents the approximation to the second derivative using standard 

centered differences, i.e., 

n rj.. ns + As,t)-2T{s,t) + T{s-As,t) 

D+-T{s,t) ^ . (2.3) 

Buttke only studied the forward case, i.e., starting from t = and, concentrating all 
the information of T at s = 0, he tried to recover the self-similar solutions we have 
described above, by imposing periodic boundary conditions at ±L. 

The scheme used by Buttke's is implicit, and he used a fixed-point iteration 
method to advance to the next time step. For this iteration to converge, a time step 
At = O(As^) is needed. Therefore, although the method is a priori unconditionally 
stable, in practice it is not. 

We have observed that an explicit finite difference scheme is equally efficient; 
more precisely, a scheme using a second-order finite difference scheme in space with 
the classical fourth order Runge-Kutta in time works well. In what follows, we have 
considered the backward case, starting at i = 1 and trying to arrive at t = 0. 

We divide [— L,+L] in N equally spaced parts, —L — sq < si < ■ ■ ■ < sjy^i < 
sn ~ +L, with 

2L 

s, = -L + iAs, i = Or--,N. 



Equation (1.9 1 is discretized as 

Tt{s,t) = T{s,t)A±D+^T{s,t), 

and the time-stepping is carried out with the classical fourth order Runge-Kutta. By 
using a fourth order Runge-Kutta, we guarantee that T o-|- T = ±1 is preserved with 
high accuracy. Nevertheless, we normalize T at every time step by doing 

ryin+l 

T"+i EE — (2 4) 

(±t::'+i o± T«+i)3 ■ 

A higher order finite difference scheme could be effortlessly implemented, although 
we have observed that this does not improve significantly the quality of our results. 

Experimentally, \At\ < 0.7As^ is found to be needed for the method to be sta- 
ble. Thus, as in Buttke's method, we have a, At = O(As^) restriction, but it is 
straightforward to advance to the next time step. 

From theorem |l.l[ we know that the self-similar solution T" has the following 
asymptotic expansion: 

T°(,s) = A^(co) - 2co^ + 0{l/.s^). (2.5) 
s 



We approximate the boundary conditions for (1.9 1 by the values of T^{s/^/t), and 
consider two types of boundary conditions: 
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1. First order boundary conditions: 



T{+L,t)=A+, 

T{-L,t)=A-. (2.6) 



2. Second order boundary conditions: 

T(+L,t)^A+-2coVf "^(+y^^ 



L 



T(-L,^) = A- + 2c„^/t^^^M), (2.7) 



normalized to be of length 1 . In ( |2.7[ ) , we use the leading order in the asymp- 
totics for 62, obtained from theorem |1.1| on page |l.l| as well. 

2.1. First order boundary conditions. We set 

T{+L,t) = A+, 

T{-L,t) = A-, yt, (2.8) 



and obtain the initial datum T(s, 1) by integrating (1.141 for s e [—L, L], using ( 1.15[ ) 



as initial condition. We have also used a fourth order Runge-Kutta, normalizing every 
T" = T(s„, 1), = ei(s„, 1) and = e^i-s,,, 1): 

rrn+l + l 
rjm+1 ^ ± gW+l = e"+-^ = — 

~ (±T"+io±T"+i)5' ^ ~ (ei+^oi 65^+1)3' ^ ~ (±e^'+io±e^+i)5 

(2.9) 

We have executed our method cq = 0.2 and different value for L. In our experiments, 
we have chosen \At\ — 0.5As^. Smaller \At\ do not improve the results, since the 
method is fourth-order accurate in time. With fixed boundary conditions, the choice 
of As is also not very important; in what follows, we show the results for L = 10 and 
L = 50, having used As — 0.01, At = —5 • 10~^. With these parameters, we have 
seen that smaller As produce virtually identical results. Nevertheless, for bigger L, 
it is convenient to use smaller As. 

To measure the quality of the results at a given i, we analyze the curvature 

In Figure [2?T| we show the curvature as a function of s and t, for L — 10 and 
L = 50. The accuracy is lost as we approach t = 0, although as one would expect, the 
bigger L is, the better the results are. However, it is remarkable the good accuracy 
with which one recovers the curvature at s = 0, even for small times, as we can see 
in Figure |2.2[ Again, this accuracy improves as we increase L. Finally, it is also 
remarkable the fact that the energy between [—L,L], 

+ L 

c^{s,t)ds, (2.10) 

-L 

calculated with the trapezoidal rule, is preserved with several precision digits (see 
Figure 2.3). This is approximately 2Lcq for all t. Since the curvature at s = is quite 



well recovered and the energy well preserved, we are able to reproduce the behavior 
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o(s.t), o -0.2. t<s[0.1,1], As-0.01, L.10 




o(s.t), o -0.2. te[0.1,1], As-0.01, L-50 




Fig. 2.1. Curvature at t € {1,0.9, ■■■ ,0.1}, with As = 0.01, L = 10,50. Although the 
results for L = 10 are very poor, for L = 50 we notice a remarkable improvement. 



c(0, t), L = 10 



c(0, t), L = 50 




O 0.1 0.2 0.3 0.4 O 0.01 0.02 0.03 0.04 



Fig. 2.2. Curvature at s = 0, with As = 0.01, L = 10, 50, together with its theoretical 
value. The accuracy improves notably when increasing L. 



Energy - 0.8. with L - 10, A s - 10 , At 5-10"^ 



O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Energy - 4, witti L - 50, As- 10 , At SIO"' 



O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Fig. 2.3. Error in the conservation of the energy, for L = 10 and L = 50, using finite 
differences and fixed extremes. Although the method is not symplectic, it preserves the energy 
very well. 



of T correctly, at least from a qualitative point of view. Indeed, in the exact solution. 
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the energy in (—00, 00) is preserved: 

c^{s,t)ds ^ 00, yt; (2.11) 



-f C30 

2/ 



this infinite energy tends to concentrate at s = 0, which causes the singularity to be 
produced. Now, we are working at [— L,i], but the finite energy is also preserved 
and it also tends to concentrate on s = 0, so we can approximate the formation of 
the singularity at i = 0. Moreover, by increasing L, we are able to recover c{0,t) for 
smaller times, improving the quality of the results. 

2.2. Second order boundary conditions. We set now 

T(+L, t)^A+- 2coVt^^^ - A+ + 2coVi ^^^''f^^''\ (2.12) 

T{-L, t)=A-+ 2coV^^^tM = A- - 2coVl^^^^^--^, (2.13) 

L L 

normalized to be of length 1 . 

Using these asymptotics as the boundary condition allows us to obtain solutions 
that arc not only qualitatively, but also quantitatively correct. The reason for this is 
that the new boundary conditions allow us to introduce energy into the system, as 
time evolves, which concentrates in order to form a cusp at the origin s — 0. 

Here, the choice of As is much more delicate. For small t and, specially, for 



big L, we have problems of aliasing, since the boundary conditions (2.121 are highly 
oscillatory. In Figure [2T4| we show the curvature at several instants in time, obtained 
using As = 0.01, As = -5 • 10~''' and L = 10 and L = 50. At t = 0.1, the curvature 
has completely degraded for L — 50. 



c(s, t), te[0.1,1], c = 0.2, L=10, As = 0.01 




Fig. 2.4. Curvature in t G {1,0.9, ■•■ ,0.1}, con As = 0.01 y L € {10,50}. The results 
tend to worsen when increasing L. 

2.3. Progressive case. In the progressive case, considered by Buttke in his 
Ph. D. thesis [6 , we start ai t ~ and go forward in time, trying to recover the 
self-similar solutions. Hence, we want to solve numerically the following initial value 
problem 

Tt(s,t) = T(s,t) A±T,,(s,t), se(-(^,+oo), 14) 
T(s,0) = Aix[o,+oo)(s) + A2x(_oo,o](s)- 
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Again, since we cannot consider the whole M, we study 



Tt{s, t) = T(s, t) A± T,,(s, i), s e [-L, 
T(s, 0) = Aix[o,+oo)(s) + A2x(_oo,o](s)- 



(2.15) 



We impose the boundary conditions T(+L,t) = T(+L,0) = and T(— L,i) = 
T(— L, 0) = A^. Dividing [— L, L] in N parts of the same length, our initial datum 
will be 



= A2, 
TO = (0,0,1), * = 
T? = Ai, 1 = 



where TO = T(s„tO) = T(s„0), 



i = 0,- 

2 
2L 



TV 



1, 



(2.16) 



N 



and T(0, 0) = (0, 0, 1), in order to be 



consistent with the backward cas e. B ecause of symmetries, T(0,t) = (0,0, 1), > 0. 

Because of part (iii) of Theorem |l.l[ (A^)3 = (A^ 
fixed, we can chose some adequate A^ and A^: 

A^ = ( \/±(l -e=F=;5'^),0,e=FT'^ 



3 = 2 Therefore, once cq is 

(2.17) 



A^ = --y/±(l -e=Fc5^),0,e 



If we executed the backward case with that cq, we would obtain those A^ and A^, 
except for a rotation around the z axis, which numerically has no relevance. 



In Figure|2!5| we have taken L = 50, 10000+1 points, i.e.. As = 0.01; M = 5-10-^ 
and Co = 0.2. We have drawn the curvature between t — 0.025 and t = 0.25, with 
increments of 0.025. It is clear that the information goes outwards with constant 
velocity. The correct theoretical curvature is shown with discontinuous strokes. We 
can appreciate that the accuracy of the results improves with time. 

c(s, t), t = 0.025i. i = 1....,10. As=10"^ 




Fig. 2.5. Curvature for different t. Progressive case with finite differences. We observe 
that the exactitude of the results improves with time; the best result is obtained at t = 0.25. 

Thus, if we want to recover the curve for a given interval at a given time, we have 
to consider a As small enough and a [—L, L] big enough. 
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If, however, we fix the boundary conditions, once the information has reached 
the boundary, it is reflected back into the domain, causing the appearance of fractal- 
like phenomena. In figure, |2.6| we have drawn the two cases of T, for cq = 0.2, 
s e [—50,50], As = 0.1, At = 5 • 10~^, t = 10. Although the approximated boundary 





-0.2 _o.5 



Fig. 2.6. Fractal formation in T 



condition is convenient for the backward case, it does not seem to be appropriate for 
the forward case, as we can see in Figure [2771 We will get back to this issue in section 



0.2, N = 1024, At = 510 , ts[1,2] 




Fig. 2.7. Progressive case, starting from t = 1, till t = 2. As we see, the approximated 
boundary condition is not adequate for the forward case. 

[3] where we consider a spectral discretization of (1.241. 

Finally, let us mention that the finite difference models have another problem in 
the forward model, which Buttke also observed. Indeed, the curvature error grows 
severely as we make At ^ (see Figure 2.8 1. In Section |3] this will also be solved. 

2.4. Evolution of X. Once T(s,t) has being calculated, it is immediate to 
recover X(s,t), since 



X,(s,t) = T(s,t) 
X(0,t) = 2cov^(0,l,0); 



(2.18) 
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c(s, t), t = 0.2, As- 0.04, c - 0.2 



A t = 2-1 



-25 -20 -15 -10 -5 O 5 10 15 20 25 



-25 -20 -15 -10 -5 O 5 10 15 20 25 



A t = 2.5-1 O • 



-25 -20 -15 -10 -5 O 5 10 15 20 25 



Fig. 2.8. Curvature noise when decreasing At 



We integrate this equation using the scheme 



As 

-j^n-l-l _ -j^n _|_ (9T" + 19T"^"'" — 5X"^^ + x"^"^) 
24 

Sn+l = Sri + As, 



(2.19) 



which, in principle, is fourth order accurate in As. Note, however, that T is only 
obtained to second order accuracy in As, resulting overall in a second order method. 
However, we use the fourth order scheme in order to reduce the accumulation of errors. 



X(s. t), 1=1,0.7,0.4,0.1,0.05,0.01, si=[-10,10] 




Fig. 2.9. Corner formation for 'K. Fixed extremes. 



We have recovered X for both boundary conditions. In Figure |2.9[ we have 
obtained T with fixed extremes and the parameters of the experiment corresponding 
to the second graphic in Figure [2?l] i.e., Euclidean case, cq = 0.2, L = 50, As = 0.01, 
Ai = -5 • 10"^, drawing only s e [-10, 10]. 

In Figure |2.11| we have drawn the X recovered from T, with the parameters 
of the experiment corresponding to Figure [2.10[ For this experiment, we had used 
the second order boundary condition, so we can see the corner formation with much 
greater clarity. Except for this, the two figures are rather similar. 
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o(s, t), ts [0.003, 0.01 5], Cj, - 0.2, L-10, As-SIO"" 

4r 



2 



-10 -5 5 10 



Fig. 2.10. Curvature at t G {0.003, 0.004, ■ ■ ■ , 0.015}, with As = 5 • 10"*. 

X(s, t), t - 1 , 0.5, 0.1 , 0.05, 0.003 




Fig. 2.11. Corner formation for Approximated boundary. 



3. Schrodinger map: Stereographic projection. We consider now the stere- 
ographic projection of T = {Ti,T2,Ts) over C: 

. = x + iy=(x,y)^(^,^). (3.1) 

We are projecting T from (0, 0, —1) into M^, identifying with C. In the EucUdean 

case, where T e there is a point on the sphere, (0.0, —1), to which no point in 
C corresponds, because the sphere is compact. Thus, we have a bijection between 
- {(0, 0, -1)} and M^. In the hyperboUc case, when T e H^, since T3 > 0, we have 
a bijection between D and H^, where 

B= {{x,y) eM.^ / x'^ <1} 

is commonly referred to as the Poincare disc. Unlike in the Euclidean, case, we do not 
need to eliminate any point from H^. The tangent vector T can be recovered from z 
by the inverse map: 

T = (Ti,T2,T'3) = ( — — — , — — — I'— — ^,7^ — o T ^2 ] • ("^•^) 
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Now, differentiating z in (3.1), together with (1.9 1 and (3.2 1, we get the following 
nonintegrable, nonlinear Schrodinger equation for z: 



2iz 



l±\z 



zl. 



2 s 



(3.3) 



The advantage of this equation, as opposed to (1.8l and (1.9l, is that the higher 
order term Zss treated implicitly, eliminating, or at least reducing significantly the 
restrictions for At. 

We can also express the other elements of the generalized Frenet-Serret formulae 
(1.11 1, i.e., c, r, ei and e.2 in terms of z. We get the generalized curvature from 



2^i 



2U. 



I±2;2±y2 l±|z|2- 

The expression for the generalized torsion is somewhat more involved: 

~ \zV± \ ^ U.J2 ■ 



(3.4) 



(3.5) 



The expression for ei follows from ei = ^T^ 



ei 



(x2 + yl)i 



1 ± a;2 ± 2/2 



1 ± a;2 ± t/2 



1 ± a;2 ± 2/2 

(3.6) 



Finally, we obtain 62 from 62 = T A± ei: 



1 



62 



(X2+ 2/2)1/2 



T2xx,y - ys{lTx'^ ±y'^) ±2xyys + Xsil ± x'^ T v'^) ±2a;ys T 2x^2/ 



1 ± a;2 -|- y2 



1 ± a;2 ± 2/^ 



1 ± a;2 ± 2/^ 
(3.7) 



3.1. Self-similar solutions. The self-similar solutions of T correspond to self- 
similar solutions of 



2iz 



l±\z 



zl 



2 s 



(3.8) 



To obtain them, note that if z is a solution of (3.8 1, so is z\{s,t) — z{Xs,X'^t), for all 
A. Taking A = t-i/2, 



z(s,t) = z{s/Vt,l) = f{s/Vt), 
where /(s) = z{s, 1). Introducing f{s/\/i) in 



Isr^/^fis/Vt) = - fis/Vi) T 2»/(5M) 1(^.(^^^))2^ 



Setting i = 1, 



r/'W-V"(.)=Fj|jfij,(/'(.))^; 
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hence 



(3.9) 



Using the same mitial conditions as before, T(0, 1) — (0, 0, 1) and Ts(0, 1) — co(l, 0, 0), 
we get 



7(o) = z(o,i) = o, 

/'(0) = z,(0,l) = 



Co 



(3.10) 



It is immediate to generahze these expressions for times other than t — 1. Defining 
g{s) = fis/Vi) = z{s,t), the EDO for g{s) is 



9"{s) = ^^/{s) ± ^ 
g{0) = /(O) = 
g'(0) = i=/'(0) 



(3.11) 



Co 

2^4 ■ 



Integrating (3.9 1, or eventually ( |3.11 1, we get the initial datum for (3.8 1. We use a 
fourth-order Runge-Kutta method to integrate (3.11 1. We obtain /(s) for s > 0. The 
values /(s) for s < are obtained by symmetry, since /(s) is antisymmetric. 

3.2. A spectral collocation method for z. We consider a semi-implicit method 



for 



Zt = IZss T 



2iz 



l±\z 



2 s ' 



(3.12) 



where we treat the linear term on the right-hand side, iZgS, implicitly, and the nonlin- 
ear term explicitly. We use a Chebyshev spectral collocation method |H], with nodes 

Si- 



Si — Lcos 



i = 0, 



and approximate z by a polynomial of the form 

N 



z{s,t) « y^^ak{t)Tk{s/L), 



fc=0 



where Tx{s) = cos(fc arccos s) is the Chebyshev polynomial of degree k. The coeffi- 
cients {flfel^o ^'"^ obtained using the fast Fourier transform (FFT) [15 . 

For the time evolution, we have chosen a second order, semi-implicit Backward 
Differentiation Formula (BDF) [22112]: 



2At 



iC/,"/^ +2AA(C/",i") -AA([/"-\r-i), (3.13) 



where we denote by Af{U,t) the nonlinear term on the right-hand side of (3.12). The 
BDF is particularly suited for this problem, as it imposes a very strong decay in the 
high frequency modes. 

Boundary conditions of Dirichlet or Neumann type can be easily implemented: 
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1. Dirichlct boundary condition: 



u{-L, t) = u,{t), J E'Loi-^r^k - Mt"+') 

k=0 "fc 



(3.14) 



2. Neumann boundary condition: 



(3.15) 



We use a spectral filter, and set equal to zero all those coefficients a^"*"^ whose modulus 



is smaller than a given e [H]. Specifically, we have made Uk = whenever ja^"^^! < 



10-^^ 

For scheme (3.131 we need two initial conditions, one at time t — , which is 
known, and the other one at t = t^. We obtain U{s,t^) « u{s,t^) using semi-implicit 
Backward Euler, and Richardson extrapolation. 

3.3. Projected second order boundary condition. In order to compare the 
finite-difference and the pseudo-spectral method, we have projected the same second 
order boundary condition we used for T: We approximate T(L,<:"'+^) as in (2.121: 

[■R„i-L^/4t"+h T"+l('r'l 

T''+\L)^A + 2coVt^^ K T"+i(L)=^ r; 

L ' (±T»+i(i)o±T"+i(i))3 



and project it over C, so 



{L) = — xT^ h i — 



.3 \L) l + Tr\L) 
A and B are computed also using the asymptotic expansion for T: 



A^T(i,io) + 2coA^^%^, 
B = {ei{L,t)-ie2{L,t))e 



where T(L,io), ei(L,to) and e2(L,io) are obtained from (3.2 1, (3.6 1 and (3.7), re- 
spectively. 

We choose the same parameters as in the finite difference case, i.e., cq = 0.2. Only 
the results for the Euclidean case will be given, since the results for the hyperbolic case 
are virtually identical for this value of cg. To measure the accuracy of our results, we 
compute the curvature as a function of space and time, and compare with the exact 
value Thus, the smaller the time that we reach with correct curvatures, the better 
we will consider the results to be. To illustrate the need for high resolution, we show 



in figure 3.1 the curvature c obtained for t = 0.04. The theoretical value is Co/^/i = 1. 
As can be seen, when N ~ 512, the accuracy is completely lost. In contrast, when 
— 2048, the error is of the order of 10~^. 
In all our numerical simulations, even for big N , we have found no restriction 
for |Ai|, so there is evidence that the method is unconditionally stable. Nevertheless, 
since the method is of order two, diminishing |Ai| can improve greatly the results. 
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-10 -5 5 10 O 1000 2000 3000 



Fig. 3.1. Curvature and spectrum of the exact z{s, 0.01), when N is too small (N — 512) 
and when N is big enough (N = 2048j. 



N = 1024, At = -10"'', 1 = 0.03,0.04,0.05 




-10 -5 O 5 10 

N = 2048, At = -10"'^, 1 = 0.03,0.04,0.05 



1 .2 h 



1.1 



0.9 




Fig. 3.2. Curvature, for |At| = 10~^. When N = 1024, the spectrum has overflowed 
and there appears noise, because of the lack of freguencies. With N — 2048, this noise have 
disappeared. 



provided N is large enough. In figure |3.2| we plot the curvature for N — 1024 and 
N = 2048, at t = 0.03, t = 0.04 y t = 0.05. For N = 1024, the noise indicates that the 
resolution is not high enough. Doubling the number of collocation points ameliorates 
the situation. Reducing the time step is not enough to remove the oscillations: using 
At = —10^®, the results at the same time instants are practically exact for N = 2048, 
but, for N = 1024 the noise has not disappeared, as we see in Figure [373] 

In order to study the cusp-formation at i = 0, we hade developed an adaptive 
methodology, both in time and space. A simple strategy is to simply duplicate the 
number of Chebyshev nodes whenever the derivative of z develops, high frequency 
components. Specifically, if we write 

N 

z,{s) = Y.^kTk{s/L), 

fe=0 

we duplicate the number of nodes when 

max |6J>2-10"''. (3.16) 
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1.3 - 
1 .2 - 



N = 1024, At = -10 , 1 = 0.03,0.04,0.05 




N = 2048, At = -10 , 1 = 0.03,0.04,0.05 



Fig. 3.3. Curvature, for \At\ = 10 Diminishing At is only useful when the number of 
frequencies is satisfactory. 



Whenever we duplicate the frequencies, we divide At by 4. 

Using this strategy, we have solved the equation starting with N = 1024, and 
Ai = -2 X 10"^ up to iV = 16384. 

Since our method is only second order in time, we would need a priori small 



\At\ in order to decrease the numerical errors. In Figure 3.4 we plot the curvature 



Adaptative method 



-4 -2 



Fig. 3.4. Adaptive method, with approximated boundary condition, from N = 1024 to 
N = 16384. 

in those instants where we have duplicated number of Chebyshev nodes, that is, at 
t = 4.91 • 10-2, t = 3.24 • 10-2, t = 1.22 • 10-^ and t = 6.09 • 10-^, and gone down up 
to t = 2.67 • 10-'^. Note that this means that the curvature has been increased by a 
factor of 20, and the energy by a factor of 400. 

3.3.1. Self-similarity boundary conditions. Since we are trying to approxi- 
mate the self-similar solutions, i.e., z{s,t) — z{s/\/t, 1), it seems natural to introduce 
this condition on the boundary. Upon differentiation, we get 21 



2t 



(3.17) 
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which can be translated into the following boundary conditions: 



2At 2t« 

z"+i(-L) - _ L 

2At ~ 21' 



z-{-L). 



This choice of the boundary condition produces good results. However, for accuracy 
reasons, for small t, we need to choose At small enough to avoid numerical artifacts 



on the boundary. This is illustrated in figure 3.5 where we show what happens on 
the boundary when At is chosen to be too large. 



N = 1024, At = -10 , te[0.2, 1] 



N = 1024, At = -5-10 , ts[0.2, 1] 



Fig. 3.5. Curvature for too big At. When At = —10 ^ , curvature explodes at s = ±L 
around t = 0.2. With At = ~5 • 10"'', there is a good improvement. 



Adaptative method 



-10 -8 



Fig. 3.6. Adaptive method, with self-similar boundary conditions, from N = 1024 till 
N = 16384. 



In Figure [XG) we have implemented an adaptive method, starting from N = 1024, 
At = — 2 • 10^^. Representing Zs{s) = J2k=o^kTk{s/L), we have duplicated the 
frequencies when 



max \bk\ > 5 • 10 



(3.18) 
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i.e., at t = 6.61 • 10"^ t = 3.23 • 10"^ t = 1.63 • 10"^ ^nd t = 8.15 • 10"'''. With 
= 16384, the results are still valid for much smaller times; in the figure we have 
also shown the curvature when t = 4.07 • 10"^, t = 3.05 • 10"^ and t = 2.67 • 10"^, 
multiplying almost by 20 the initial curvature. Observe that, in comparison with 
Figure |3.4[ the results of Figure 3^ are much cleaner. 

3.4. Radiation boundary condition. We derive another boundary condition, 
that aims at capturing the correct flow of energy through the boundary at s = ±L. 
For this, we consider the Hasimoto transform, ip = cexp (i r(s',t)ds'). From (3.4 1 
and (3.5 1, it follows that 



2U, 



l± z 



■ exp 



\±l + x'^ + y'^ 



xj + yj 



ds' 



Now, 



Vs^s 



d f ys 

= — arctan — 

OS \ Xo 



so 



2U., 



1± |z|2 

2\z,\ 
1± |z, 

l±\z\ 



exp 
2 exp 
■ exp 



{■r 



2{yxs - xys) ^ d_ 
±1 + + ds 



i arctan 



2{yxs - xys) 



±1 + x2 + j/2 



i arctan 
exp 



arctan | — 

x,s{0) 

i arctan 



ds' 



exp 



2{yxs - xys) 



±l + a;2 + y2 



Xs{0) 



(3.19) 



Co 



At s = ±L, we know that ■0(±L,i) = • From this, it follows that 



Co 



resulting 



2zs 



l±\z\ 



exp 



2{yxs " xys) 
±l + x'^ + y2 



exp 



-i arctan 



Xs{0) 



s=±L 



Zs{±L,t) 



2 Vt 



exp 



2{yxs - xys) 
±l + a;2 + y2 ^ 



exp 



i arctan 



Xs(0) 



s=±L 



We can give the boundary condition as 
"l±kP CO 



zT\±L) 



exp 



2{yxs - xys) 
±l + a;2 + y2 * 



exp 



i arctan 



(3.20) 

x,(0) 



t=t" 



which makes the scheme be first order. To obtain a second order scheme, we do 
12 , r ou.^ ^ 1 r /y^(o) 



(±L) = 2 



1± Co ,Z 

e 1' exp 



2 ^/t 
1±NP CO 



2{yxs ~ xys) , 
±l + a;2 + y2 



exp 



i arctan 



e 4* exp 



2{yxs -~ xys) , 

— ^ T^ds exp 

±l + a::2 + ?;2 



i arctan 



Xs{0) 

VsjO) 

Xs{0) 



(3.21) 



s = ±L 



s=±L 



A numerical study of the self-similar solutions of the Schrodinger Map 



23 



Notice that for unperturbed solutions or for solutions perturbed with an even pertur- 
bation, we have 



exp 



i arctan 



XsiO) 



We have found the scheme to be unconditionally stable, and that the results in the 
backward case are as good as those obtained with the boundary condition described 
m section [Ol] What is more important, this boundary condition seems to work 
extraordinarily well for the progressive case. In Figure 3.7 considering t'^ — 1, we 
have solved the equation up to < = 2. The results in this case are excellent. 



. 0.2, N = 1024, A t = 1 10" 



ts [1 , 2] 



Fig. 3.7. Forward case, using the third boundary condttions. In comparison with figure 

\2. 7[ the results are excellent. 



3.5. Forward case in time. The pseudo-spectral method, together with the 
last boundary condition, allows to give a complete treatment of the forward case, 
unlike with the finite differences. 

We consider the following initial value problem: 

2iz „ 



^z{s,t) = aisx[o,+oo)(s) + a^sX(-oo,o](s). 

Once Co is fixed, we can consider take and a? as the respective projections of A^(co) 
and A^(co) used in (2.171, i.e.. 



±(1 - eT^o'^) 



a2 = - 



±(1 - eT'^o^) 



(3.22) 



Therefore, discretizing [— L, L] in = Lcos{iTT/N), with i = 0, - ■ ■ ,N, the numerical 
initial datum is 



N 
' 2 



1}, 



'zO = z(s„0) EE a+, ie{0, 
4/2 = z(0,0) = 0, 
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This boundary condition is not adequate for large t; when the information reaches 
the boundary, it is reflected back into the domain, creating a fractal phenomenon, 
as illustrated in figure |3.8| where we show the case cq = 0.2, N = 16384, L = 50 
and t = 10. If, however, we consider for the same parameters a much smaller time. 




1.15 I ' ' ' ' ' ' 

-0.3 -0.2 -0.1 O 0.1 0.2 0.3 



Fig. 3.8. Fractal creation 



for instance t = 0.3, we observe in Figure |3.9| that the information, even if it has 
already reached the boundary, has only partially rebounded. Thus, in spite of the 
great noise at the extremes, if we make a zoom of the central subinterval [—10,10], 
the achieved accuracy is notorious. Therefore, we will consider the portion of the 



ss [-50, 50] 



ss[-10, 10] 





Fig. 3.9. Starting from t = 0, with At —10 ^ , att = 0.3 the information has reached the 
boundary and partially rebounded. Nevertheless, m s £ [—10, 10], the error of the curvature 
is smaller that 9 • 10^*, so that portion of curve is well suited to be our new initial datum. 

curve corresponding to that subinterval as our new initial datum. 

Now, we know z(s, 0.3) at Si — 50 cos(i7r/16384), j = 0, • • • , 16384, and have to 
interpolate spectrally z(s,0.3) in the new initial nodes Si, which we have chosen as 



Si — 10 cos 



1024 



i = 0, • • • , 1024. 



Therefore, we have to interpolate z in 1024+ 1 points belonging to [—10, 10]; this can 



A numerical study of the self-similar solutions of the Schrodinger Map 25 

be done directly from the expression for z as a function of the coefficients a^: 




Fig. 3.10. Taking as initial datum the right side of figure \ 3.9[ we have advanced till 
t — 1.5, with At — 10~* and the third boundary condition. At the right side, we have drawn 
the curve obtained at t — 1.5. 

In Figure |3.10| we have calculated the evolution of z with that new initial datum 
and the third boundary condition in s = ±10. Although there appears some noise, 
because the new initial datum is not exact, the results are acceptable. 

To finish, let us mention also that, unlike in the finite difference case, there are 
no problems when At ^ in any of the two stages. 

4. A stability result. In the previous sections, we have described several schemes 
to approximate numerically the self-similar solutions of 

Xt - X, A± X,, (4.1) 

and 

Ti = TA±T,„ (4.2) 

which are characterized by 

c(.,t)=|, r{s,t) = ^^. (4.3) 

Since we could not work with all M, we have bounded ourselves to s G [— L, L], 
considering different boundary conditions at s = ±L. It is clear that these boundary 
conditions, from a numerical point of view, will never be completely exact; even in 
the best of the cases there will be some errors, although they may be as small as the 
machine precision. Thus, we can consider all our experiments to be perturbation of 
the exact solutions for X, T and the projection of T, z. 

In the previous sections, we have seen that the stability of all the experiments is 
good, no matter wether we consider very rough boundary conditions, as fixing T(s, t) 
at ±L or more elaborated ones. Thus, before concluding this paper, it is interesting 
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to mention for completeness some recent results by Banica and Vega [3l H], that 
guarantee theoretically this stability, at least for small perturbations. 

Equations (4.1 1 and (4.2 1 are closely related to the cubic Schrodinger equation or 
non-linear Schrodinger equation (NLS). Indeed, by means of Hasimoto's transforms 

ip{s,t) = c{s,t)exp (i [ T{s',t)ds'\ (4.4) 



we obtain the NLS [2D] [B], 



(4.5) 



The term A{t) can be immediately absorbed by means of a change of variable "if 



i?/'t + V'ss±-|V'PV==0. 



(4.6) 



The case with the + sign is known as the focussing case, the cubic Schrodinger equa- 
tion being denote NLS^; this corresponds to the Euclidean case of (1.8 1 and (1.9 1. 
With the — sign, corresponding to the hyperbolic case of (1.8) and (1.9), we have 
the defocussing case, the Schrodinger cubic being denoted as NLS^. The Schrodinger 
cubic equation appears in many contexts |12j . as certain non- linear optics phenomena, 
wave packets in water and plasma, etc... 



Coming back to (4.5), in our current problem, from (4.3), 



_£0_ ,sV4t 



(4.7) 



If we choose A{t) = —'f, we get a solution of (4.5), with ?A(0,s) — V^cqS, being 6 
Dirac's distribution. Thus, going from X and z to means to trivialize the problem 
in a certain way, since, for every t > 0, we know the explicit equation. 

Now, the solutions (4.7) of the Schrodinger cubic equation (4.5) that we are 
considering satisfy trivially 



Therefore, from this point of view, our solutions have infinite energy. Nevertheless, 
Banica and Vega [2111] have proved that, under certain renormalizations, the solutions 
have finite energy. Indeed, starting from 



we make the change u{s, t) = ip{\/2s, 2t) in order to absorb the constant i, obtaining 



iut{s,t) + Uss{s,t) ± \u{s,t)\ 



2t 



u{s,t) = 0. 



(4.8) 



Applying to this last expression the following conformal transform 
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and evaluating in (s, t) — (|, 1), we get 



vt{s,t) = ~ivss{s,t) =F - 
Therefore, when < t < to, u is solution of 



iut{s,t) + u,s{s,t) ± {\u{s,t)\^ ~ ^ju{s,t) = 

u{ta,s) : 

if and only if w is a solution, when l/to < t < oo, of 
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(4.10) 



0, 



vt{s,t) = -iVss{s,t)T -\\vis,t)\ 



Co 



2 eg 



(4.11) 



w(s,t), 



(4.12) 



with wo(s) = r^^ui(s). Notice that, due to the conformal transform, the case i ^ 0+ 
gets transformed into t +oo. Moreover, our exact unperturbed solutions are now 
constant solutions of u. 

There is an energy naturally associated to (4.121, 

J \vs{s,t)\'dsTl^ J (Ks,i)p-f) ds. (4.13) 



Therefore, if w is a solution of (4.121, we get 



d_ 
di 



J (\v{s,t)\' 



ds = 0. 



(4.14) 



This last equation implies, in the defocussing situation corresponding to the hyperbolic 
case, that the energy does not grow when t ^ oo. As a consequence of (4.14), Banica 
and Vega proved in [3] the following theorem for the defocussing case: 
Theorem 4.1. 

For every to > and for every vq e Ti.^ , there exists a unique solution of the 
initial value problem {4-- 12), with 

v-^ eC{{i/to,^),n'). 

Banica and Vega also proved that 

1 



lim inf 

t — ^oo t 



Ks,t)\'-'f 



ds = 0, 



which implies that u, in the defocussing case, when < t < to, satisfies 

limjnf\\t\u{t)f-clh^O. 

This last expression is already a stability result for the singular solution 

of (liTTl). 



/Wo 



-e 4 
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Finally, in the defocussing case, writing in terms of the geometric quantities c and 



T the corresponding energy to (4.4 1, which is a solution of (4.5 1, we get 

Ids 



and, hence, 



f2 r+co / ; 

/ + 00 
[tc^{s,t) - cl]'^ds, (4.15) 
OO 



16^2 



d 



1 f^°° 

m - J ^ - col'rfs = (4.16) 



and 



liminf|lt|c|2-cg|l2 = 0. (4.17) 

Observe that our solutions are precisely such that E{t) = 0, for all t > 0. 

In the focussing setting (i.e. the case of the sphere) the stability is much more 
delicate. In [?^, and under a smallness assumption in the curvature cq, they construct 



a global solution such that E{t) given in (4.131 is finite. The main difficulty comes 



from the long range character of the non- linear potential that appears in (4.121. This 
implies the existence of a logarithmic phase that has as a consequence the non ex- 
istence of the limit at infinity for the solutions. However they also prove, this time 
in the focussing case, that this logarithmic divergence disappears when the tangent 
vector T is computed, so that there is stability for T. Remember that in order to 
compute T one has to integrate once the curvature and twice the torsion. Therefore 
thanks to the oscillations the integrals converge without any difficulty. 

5. Conclusions. In this paper, we have tried to reproduce numerically the be- 
havior of the self-similar solutions of 

Xt = X, A± X,, (5.1) 

and 

Ti = TA±T,„ (5.2) 
which develop a singularity at finite time. These solutions, characterized by 

c(s,i)=y, 

form a one-parameter family, where cq is precisely the family parameter. 



The singularities happen at t = 0, going backwards in time, but, since both (5.1 1 



and (5.2 1 are time-reversible, we could consider an equivalent problem, where we 
advance in time and the singularity happens. 

In Section [2] we have given a finite-difference scheme to study the self-similar 
solutions of T. T determines completely X, except for a constant, determined by 
X(0,t) = 2cov^(0, 1, 0). We have studied mainly the backward case: starting from 
t = 1, we have tried to reproduce the formation of the singularity. Since we cannot 
consider all M, we have bounded ourselves to s G [—L.L], being necessary to give 
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boundary conditions a,t s — ±L. From the asymptotics of Theorem |1.1| for T(s, 1), 
we have deduced two approximated boundary conditions, depending on whether we 
choose the leading term of T(s, 1), or the first two terms in the expansion. 

Considering T(s, 1) constant at s = ±L gives us good results from a qualitative 
point of view: the energy is preserved with several precision digits and we get also 
the approximated value of c(0, t), even for small times. Thus, the finite energy tends 
to concentrate on the origin, approximating the formation of the singularity ai t — 0. 
In the exact solution, the energy was infinite for all t, but the behavior is the same: 
all the energy tends to concentrate at s = 0. 

Since the bigger is L, the smaller are the times for which we recover c(0, t), there 
seems to be evidence that we could approximate the exact problem by making L tend 
to infinity. It would be very interesting to prove this analytically. 

The second boundary condition is obtained considering the first non-constant 
term in the asymptotics of T(s, 1), and that T(s,t) — T(s/\/t, 1). For not too big L 
and provided that As is small enough, it allows us to recover the solutions with big 
accuracy even for small times; moreover, not only for s = 0, but for all s. Nevertheless, 
since \At\ = 0{As^), this is quite expensive from a computational point of view. 

When L is bigger, a uniform distribution of the nodes is not adequate, because 
there is a lack of resolution for s near the boundary. That suggests that we use an 
alternative node distribution, given by the Chebyshev nodes, which are distributed 
in a much more suitable way. Since, with an explicit scheme, we would have now a 
\At\ = 0{N~'^) restriction, we have projected stereographically T over C, obtaining 



We have implemented an implicit-explicit method in time, with a pseudo-spectral 
method in space, considering different boundary conditions. This method has several 
interesting advantages: its stability is much better for a given TV; it allows consid- 
ering much bigger L; it allows implementing boundary conditions that do not need 
external information and X can be recovered immediately with spectral accuracy. 
Since designing an adaptive version, i.e., refining the grid when necessary, is natural 
and straight-forward, the method is globally much more efficient. Lastly, it allows to 
make a full treatment of the forward case: starting from a singular datum at t = 0, 
we can recover the solutions in two stages, by adapting at i = e adequate boundary 
conditions. 

Finally, let us mention that all the experiments we have done give evidence of the 



stability of equations (5.11, (5.2 1 and (5.4 1 from a numerical point of view. Indeed 



when giving the boundary conditions at s = ±L, it is impossible to do it exactly 
and we will always be introducing some perturbations; when, for instance, we fix 
T(s,t) to be constant at s = ±L, these perturbations will be quite big. Even if we 
calculated the exact value of X, T or ^ at s = ±L, there would always be a tiny 
error attributable to machine precision. Therefore, in Section [4] we have mention a 
recent stability result by Banica and Vega, concerning the self-similar solutions we 
have studied. This offers, somehow, a theoretical support to this paper. 
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